A Flexible Bayesian Model for Studying Gene–Environment Interaction

An important follow-up step after genetic markers are found to be associated with a disease outcome is a more detailed analysis investigating how the implicated gene or chromosomal region and an established environment risk factor interact to influence the disease risk. The standard approach to this study of gene–environment interaction considers one genetic marker at a time and therefore could misrepresent and underestimate the genetic contribution to the joint effect when one or more functional loci, some of which might not be genotyped, exist in the region and interact with the environment risk factor in a complex way. We develop a more global approach based on a Bayesian model that uses a latent genetic profile variable to capture all of the genetic variation in the entire targeted region and allows the environment effect to vary across different genetic profile categories. We also propose a resampling-based test derived from the developed Bayesian model for the detection of gene–environment interaction. Using data collected in the Environment and Genetics in Lung Cancer Etiology (EAGLE) study, we apply the Bayesian model to evaluate the joint effect of smoking intensity and genetic variants in the 15q25.1 region, which contains a cluster of nicotinic acetylcholine receptor genes and has been shown to be associated with both lung cancer and smoking behavior. We find evidence for gene–environment interaction (P-value = 0.016), with the smoking effect appearing to be stronger in subjects with a genetic profile associated with a higher lung cancer risk; the conventional test of gene–environment interaction based on the single-marker approach is far from significant

A regression model for risk difference estimation in population-based case-control studies clarifies gender differences in lung cancer risk of smokers and never smokers

BACKGROUND: Additive risk models are necessary for understanding the joint effects of exposures on individual and population disease risk. Yet technical challenges have limited the consideration of additive risk models in case-control studies. METHODS: Using a flexible risk regression model that allows additive and multiplicative components to estimate absolute risks and risk differences, we report a new analysis of data from the population-based case-control Environment And Genetics in Lung cancer Etiology study, conducted in Northern Italy between 2002-2005. The analysis provides estimates of the gender-specific absolute risk (cumulative risk) for non-smoking- and smoking-associated lung cancer, adjusted for demographic, occupational, and smoking history variables. RESULTS: In the multiple-variable lexpit regression, the adjusted 3-year absolute risk of lung cancer in never smokers was 4.6 per 100,000 persons higher in women than men. However, the absolute increase in 3-year risk of lung cancer for every 10 additional pack-years smoked was less for women than men, 13.6 versus 52.9 per 100,000 persons. CONCLUSIONS: In a Northern Italian population, the absolute risk of lung cancer among never smokers is higher in women than men but among smokers is lower in women than men. Lexpit regression is a novel approach to additive-multiplicative risk modeling that can contribute to clearer interpretation of population-based case-control studies

Risk of miscarriage with bivalent vaccine against human papillomavirus (HPV) types 16 and 18: pooled analysis of two randomised controlled trials

Objective To assess whether vaccination against human papillomavirus (HPV) increases the risk of miscarriage

Lung Cancer and Occupation in a Population-based Case-Control Study

The authors examined the relation between occupation and lung cancer in the large, population-based Environment And Genetics in Lung cancer Etiology (EAGLE) case-control study. In 2002–2005 in the Lombardy region of northern Italy, 2,100 incident lung cancer cases and 2,120 randomly selected population controls were enrolled. Lifetime occupational histories (industry and job title) were coded by using standard international classifications and were translated into occupations known (list A) or suspected (list B) to be associated with lung cancer. Smoking-adjusted odds ratios and 95% confidence intervals were calculated with logistic regression. For men, an increased risk was found for list A (177 exposed cases and 100 controls; odds ratio = 1.74, 95% confidence interval: 1.27, 2.38) and most occupations therein. No overall excess was found for list B with the exception of filling station attendants and bus and truck drivers (men) and launderers and dry cleaners (women). The authors estimated that 4.9% (95% confidence interval: 2.0, 7.8) of lung cancers in men were attributable to occupation. Among those in other occupations, risk excesses were found for metal workers, barbers and hairdressers, and other motor vehicle drivers. These results indicate that past exposure to occupational carcinogens remains an important determinant of lung cancer occurrence

A Case-Control Study of Peripheral Blood Mitochondrial DNA Copy Number and Risk of Renal Cell Carcinoma

Background: Low mitochondrial DNA (mtDNA) copy number is a common feature of renal cell carcinoma (RCC), and may influence tumor development. Results: from a recent case-control study suggest that low mtDNA copy number in peripheral blood may be a marker for increased RCC risk. In an attempt to replicate that finding, we measured mtDNA copy number in peripheral blood DNA from a U.S. population-based case-control study of RCC. Methodology/Principal Findings: Relative mtDNA copy number was measured in triplicate by a quantitative real-time PCR assay using DNA extracted from peripheral whole blood. Cases (n = 603) had significantly lower mtDNA copy number than controls (n = 603; medians 0.85, 0.91 respectively; P = 0.0001). In multiple logistic regression analyses, the lowest quartile of mtDNA copy number was associated with a 60% increase in RCC risk relative to the highest quartile (OR = 1.6, 95% CI = 1.1–2.2; Ptrend = 0.009). This association remained in analyses restricted to cases treated by surgery alone (OR Q1 = 1.4, 95% CI = 1.0–2.1) and to localized tumors (2.0, 1.3–2.8). Conclusions/Significance: Our findings from this investigation, to our knowledge the largest of its kind, offer important confirmatory evidence that low mtDNA copy number is associated with increased RCC risk. Additional research is needed to assess whether the association is replicable in prospective studies

Short term persistence of human papillomavirus and risk of cervical precancer and cancer: population based cohort study

Objective To evaluate the cumulative incidence of cervical intraepithelial neoplasia II or worse (grade II+) or cervical intraepithelial neoplasia grade III+ after short term persistence of prevalently detected carcinogenic human papillomavirus (HPV)

Population Substructure and Control Selection in Genome-Wide Association Studies

Determination of the relevance of both demanding classical epidemiologic criteria for control selection and robust handling of population stratification (PS) represents a major challenge in the design and analysis of genome-wide association studies (GWAS). Empirical data from two GWAS in European Americans of the Cancer Genetic Markers of Susceptibility (CGEMS) project were used to evaluate the impact of PS in studies with different control selection strategies. In each of the two original case-control studies nested in corresponding prospective cohorts, a minor confounding effect due to PS (inflation factor λ of 1.025 and 1.005) was observed. In contrast, when the control groups were exchanged to mimic a cost-effective but theoretically less desirable control selection strategy, the confounding effects were larger (λ of 1.090 and 1.062). A panel of 12,898 autosomal SNPs common to both the Illumina and Affymetrix commercial platforms and with low local background linkage disequilibrium (pair-wise r2<0.004) was selected to infer population substructure with principal component analysis. A novel permutation procedure was developed for the correction of PS that identified a smaller set of principal components and achieved a better control of type I error (to λ of 1.032 and 1.006, respectively) than currently used methods. The overlap between sets of SNPs in the bottom 5% of p-values based on the new test and the test without PS correction was about 80%, with the majority of discordant SNPs having both ranks close to the threshold. Thus, for the CGEMS GWAS of prostate and breast cancer conducted in European Americans, PS does not appear to be a major problem in well-designed studies. A study using suboptimal controls can have acceptable type I error when an effective strategy for the correction of PS is employed

Impact of human papillomavirus (HPV) 16 and 18 vaccination on prevalent infections and rates of cervical lesions after excisional treatment

BackgroundHuman papillomavirus vaccines prevent human papillomavirus infection and cervical precancers. The impact of vaccinating women with a current infection or after treatment for an human papillomavirus-associated lesion is not fully understood.ObjectivesTo determine whether human papillomavirus-16/18 vaccination influences the outcome of infections present at vaccination and the rate of infection and disease after treatment of lesions.Study DesignWe included 1711 women (18−25 years) with carcinogenic human papillomavirus infection and 311 women of similar age who underwent treatment for cervical precancer and who participated in a community-based trial of the AS04-adjuvanted human papillomavirus-16/18 virus-like particle vaccine. Participants were randomized (human papillomavirus or hepatitis A vaccine) and offered 3 vaccinations over 6 months. Follow-up included annual visits (more frequently if clinically indicated), referral to colposcopy of high-grade and persistent low-grade lesions, treatment by loop electrosurgical excisional procedure when clinically indicated, and cytologic and virologic follow-up after treatment. Among women with human papillomavirus infection at the time of vaccination, we considered type-specific viral clearance, and development of cytologic (squamous intraepithelial lesions) and histologic (cervical intraepithelial neoplasia) lesions. Among treated women, we considered single-time and persistent human papillomavirus infection, squamous intraepithelial lesions, and cervical intraepithelial neoplasia 2 or greater. Outcomes associated with infections absent before treatment also were evaluated. Infection-level analyses were performed and vaccine efficacy estimated.ResultsMedian follow-up was 56.7 months (women with human papillomavirus infection) and 27.3 months (treated women). There was no evidence of vaccine efficacy to increase clearance of human papillomavirus infections or decrease incidence of cytologic/histologic abnormalities associated with human papillomavirus types present at enrollment. Vaccine efficacy for human papillomavirus 16/18 clearance and against human papillomavirus 16/18 progression from infection to cervical intraepithelial neoplasia 2 or greater were −5.4% (95% confidence interval −19,10) and 0.3% (95% confidence interval −69,41), respectively. Among treated women, 34.1% had oncogenic infection and 1.6% had cervical intraepithelial neoplasia 2 or greater detected after treatment, respectively, and of these 69.8% and 20.0% were the result of new infections. We observed no significant effect of vaccination on rates of infection/lesions after treatment. Vaccine efficacy estimates for human papillomavirus 16/18 associated persistent infection and cervical intraepithelial neoplasia 2 or greater after treatment were 34.7% (95% confidence interval −131, 82) and −211% (95% confidence interval −2901, 68), respectively. We observed evidence for a partial and nonsignificant protective effect of vaccination against new infections absent before treatment. For incident human papillomavirus 16/18, human papillomavirus 31/33/45, and oncogenic human papillomavirus infections post-treatment, vaccine efficacy estimates were 57.9% (95% confidence interval −43, 88), 72.9% (95% confidence interval 29, 90), and 36.7% (95% confidence interval 1.5, 59), respectively.ConclusionWe find no evidence for a vaccine effect on the fate of detectable human papillomavirus infections. We show that vaccination does not protect against infections/lesions after treatment. Evaluation of vaccine protection against new infections after treatment and resultant lesions warrants further consideration in future studies

Associations of CDH1 germline variant location and cancer phenotype in families with hereditary diffuse gastric cancer (HDGC)

INTRODUCTION: Hereditary diffuse gastric cancer (HDGC) is a cancer syndrome associated with variants in E-cadherin (CDH1), diffuse gastric cancer and lobular breast cancer. There is considerable heterogeneity in its clinical manifestations. This study aimed to determine associations between CDH1 germline variant status and clinical phenotypes of HDGC. METHODS: One hundred and fifty-two HDGC families, including six previously unreported families, were identified. CDH1 gene-specific guidelines released by the Clinical Genome Resource (ClinGen) CDH1 Variant Curation Expert Panel were applied for pathogenicity classification of truncating, missense and splice site CDH1 germline variants. We evaluated ORs between location of truncating variants of CDH1 and incidence of colorectal cancer, breast cancer and cancer at young age (gastric cancer at \u3c40 or breast cancer \u3c50 years of age). RESULTS: Frequency of truncating germline CDH1 variants varied across functional domains of the E-cadherin receptor gene and was highest in linker (0.05785 counts/base pair; p=0.0111) and PRE regions (0.10000; p=0.0059). Families with truncating CDH1 germline variants located in the PRE-PRO region were six times more likely to have family members affected by colorectal cancer (OR 6.20, 95% CI 1.79 to 21.48; p=0.004) compared with germline variants in other regions. Variants in the intracellular E-cadherin region were protective for cancer at young age (OR 0.2, 95% CI 0.06 to 0.64; p=0.0071) and in the linker regions for breast cancer (OR 0.35, 95% CI 0.12 to 0.99; p=0.0493). Different CDH1 genotypes were associated with different intracellular signalling activation levels including different p-ERK, p-mTOR and β-catenin levels in early submucosal T1a lesions of HDGC families with different CDH1 variants. CONCLUSION: Type and location of CDH1 germline variants may help to identify families at increased risk for concomitant cancers that might benefit from individualised surveillance and intervention strategies

Phase I Metabolic Genes and Risk of Lung Cancer: Multiple Polymorphisms and mRNA Expression

Polymorphisms in genes coding for enzymes that activate tobacco lung carcinogens may generate inter-individual differences in lung cancer risk. Previous studies had limited sample sizes, poor exposure characterization, and a few single nucleotide polymorphisms (SNPs) tested in candidate genes. We analyzed 25 SNPs (some previously untested) in 2101 primary lung cancer cases and 2120 population controls from the Environment And Genetics in Lung cancer Etiology (EAGLE) study from six phase I metabolic genes, including cytochrome P450s, microsomal epoxide hydrolase, and myeloperoxidase. We evaluated the main genotype effects and genotype-smoking interactions in lung cancer risk overall and in the major histology subtypes. We tested the combined effect of multiple SNPs on lung cancer risk and on gene expression. Findings were prioritized based on significance thresholds and consistency across different analyses, and accounted for multiple testing and prior knowledge. Two haplotypes in EPHX1 were significantly associated with lung cancer risk in the overall population. In addition, CYP1B1 and CYP2A6 polymorphisms were inversely associated with adenocarcinoma and squamous cell carcinoma risk, respectively. Moreover, the association between CYP1A1 rs2606345 genotype and lung cancer was significantly modified by intensity of cigarette smoking, suggesting an underling dose-response mechanism. Finally, increasing number of variants at CYP1A1/A2 genes revealed significant protection in never smokers and risk in ever smokers. Results were supported by differential gene expression in non-tumor lung tissue samples with down-regulation of CYP1A1 in never smokers and up-regulation in smokers from CYP1A1/A2 SNPs. The significant haplotype associations emphasize that the effect of multiple SNPs may be important despite null single SNP-associations, and warrants consideration in genome-wide association studies (GWAS). Our findings emphasize the necessity of post-GWAS fine mapping and SNP functional assessment to further elucidate cancer risk associations